
img = imread('1.bmp');
img = img(:,:,1);
I = img(:);
[m,n] = size(img);
L = 22;
matrix = getSampMatrix(n,L);
A = double(matrix);
   
b = A * double(I);

maxIter = 200*m ;
%x01 = A\b;
x0 = zeros(m*n,1);
tol = 10;
[x res] = reconART(b,A,tol,x0,maxIter);
